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Dear Editors, 

Middle East Respiratory Syndrome (MERS) has spread to 
more than 23 countries since the first case of a MERS 
coronavirus (MERS CoV) infection was confirmed in Sep 
2012 [1], Because MERS CoV has led to mortality as high 
as 40%, it has become one of the most important emerging 
pathogens threatening human health. Recently, the outbreak 
of MERS CoV, which resulted in 172 infections and 27 
deaths in South Korea, was the largest outside Saudi Arabia. 
A 44-year-old man, had been exposed to the first confirmed 
MERS case who shared a ward with the man’s father in this 
outbreak, travelled from South Korea to China, becoming 
the first imported MERS case reported in China. Here, we 
report the sequencing and analysis of the MERS CoV ge¬ 
nome from this case. Viral genome sequencing provides 
information vital to tracing the pathogenic source, predict¬ 
ing viral virulence, and following the pattern of viral evolu¬ 
tion, variation, and transmission [2-4], 

Twenty-four primer pairs, specific to conserved regions, 
were used to amplify the complete genome of MERS CoV 
by PCR. The resulting PCR products had an average length 
of 1,500 base pairs (bp). The PCR products were recovered 
and sequenced by the Sanger dideoxy method. The genomic 
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sequence was assembled with DNAstar 7.0 software and 
deposited in GenBank with accession no. KT036372. The 
deposited sequence was 29,928 bp, which was -180 bp 
shorter than the complete genomic sequence, as the tech¬ 
nique of 5'- and 3'-RACE (rapid-amplification of cDNA 
ends), which amplifies genomic ends, was not applied. 

The sequence we obtained: KT036372, differed slightly 
from other MERS CoV sequences submitted. It differed in 
four nucleotides from the updated sequence of KT006149.2 
provided by China Centers for Disease Control (CDC) from 
the same patient; and in 10 nucleotides from sequence 
KT029139 submitted by Korea Centers for Disease Control 
and Prevention. The two sequences used for comparison, 
KT006149.2 and KT029139, differed from each other by 12 
nucleotides. Our analysis shows that the KT036372 se¬ 
quence is reliable, sharing greater similarity to the 
KT029139 sequence from South Korea. The differences 
among the three sequences (Table SI) might have resulted 
from: viral evolution and mutation; the collection date 
(KT036372 on May 28; KT006149 on May 27); sample 
origin (KT036372 and KT006149 from clinical sample; 
KT029139 from virus culture); and the sequencing method 
(KT036372 by Sanger dideoxy method; KT006149 by Ion 
Torrent platform and Sanger dideoxy method; and 
KT029139 using the Illumina platform). 

Complete alignment of the genomic sequence in Gen- 
Bank indicated that KT036372 shared a high similarity of 
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Figure 1 Phylogenetic tree of MERS CoV sequence generated by the Neighbor Joining method using MEGA6.0 and tested with 1,000 bootstrap replicates, 
a, Phylogenetic tree of MERS CoV based on genomic sequence, b. Phylogenetic tree of MERS CoV based on ORFla gene sequence. H, The MERS isolate 
determined in this study; A, The MERS strain from the same patient sequenced by China CDC; ▼, The sequence submitted by South Korea. Phylogentic 
trees show that MERS-CoV sequences identified in Korea and China cluster closest with viruses obtained from a recent outbreak in Riyadh in 2015 (KSAa). 
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99.53%-99.92% with the strains isolated in recent years, the 
highest identity being with the Saudi Arabian isolate se¬ 
quenced in Mar 2015. The S gene of the virus, which is as¬ 
sociated with viral factors for host recognition, also shared 
high similarity of 98.97%-99.93% with the S gene of the 
MERS CoV sequences stored in the database, with highest 
similarity to KT026453. 

We found three novel variants, with single amino acid 
substitutions: D977G, T1833I, and G6896S, which differ 
from the strains reported in the database. These mutations 
are located in the poly-protein lab (ORFlab) that is cleaved 
into 16 nonstructural proteins (NSP1-16). D977G and 
T1833I are situated in the nonstructural protein 3 (NSP3), 
predicted to be a protease responsible for intracellular post- 
translational processing of viral proteins; while the mutation 
G6896S is located in the nonstructural protein 16 (NSP16), 
whose function is unknown. However, no new mutation has 
emerged in the proteins closely relating to virus assembly 
and virulence, such as the S, E, M, and N proteins. 

Phylogenetic analysis of the complete genomic sequence 
further demonstrated that KT036372 as well as KT029139, 
belonged to clade B (Figure la), like the KSAa strains re¬ 
cently isolated from Saudi Arabia. Sufficient sequence di¬ 
vergence was not seen for them to be assigned to a new 
clade. Remarkably, there were two different secondary 
branches KSAa and KSAb, suggesting that two different 
origins of MERS CoV may prevail in Saudi Arabia. The 
phylogenetic tree of ORFla gene (Figure lb) was closer to 
that of the full-length genome. Nevertheless, S gene phylo¬ 
genetic tree was similar to that of the genome, but distinct in 
that KSAa and KSAb were clustered together (Figure SI). 

The full-length genomic sequence exhibited high similar¬ 
ity above 95% with MERS CoV sequences found in the 
database. No recombination events were identified in the 
current epidemic strains by either the Simplot or the RDP 
software methods. 

In conclusion, the sequence of the MERS CoV isolate 
from the first imported case in China’s Guangdong prov¬ 
ince, was highly similar to that of recent isolates from the 
Saudi Arabia. We identified three new mutations in the 
polyprotein ORFlab; however, no variation in genes asso¬ 
ciated with viral assembly or the virulence-linked S, E, M 


and N genes were observed, nor was there any apparent 
recombination in the genome. These results further suggest 
that the virulence of the MERS CoV did not increase. Based 
on current reports, mortality was -10%, which is far lesser 
than the 40% previously reported. Simultaneously, the epi¬ 
demiological characteristics of the virus suggest that it is 
still limited to human-to-human transmission, with no sign 
of a large-scale outbreak. 

As the case reported here illustrates, globalization and 
frequent international travel can pose a threat to human 
health worldwide, and exposure to the infectious agent can 
occur unpredictably. Hence, excellent technical preparations 
and mandatory close monitoring are critical for rapid detec¬ 
tion of pathogenic viruses like MERS CoV, and to reduce 
the risk of a pandemic. We were able to achieve this by 
identification and analysis of the viral genome from the 
only imported case of MERS CoV, with no contacts being 
infected in Guangdong up to now [5]. 
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Figure SI Phylogenetic tree of S gene of MERS CoV built with Neighbor Joining method using MEGA6.0 and tested with 1,000 bootstrap replicates. 50 
strains with high similarity isolated from Saudi Arabia and United Arab Emirates in 2013-2014 were compressed in the phylogenetic tree of S gene. 
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